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We investigate the 16-band d-p model for iron pnictide superconductors in the presence of 
the electron-phonon coupling g with the orthorhombic mode which is crucial for reproducing 
the recently observed ultrasonic softening. Within the RPA, we obtain the ferro-orbital order 
below Tq which induces the tetragonal-orthorhombic structural transition at T s = Tq, together 
with the stripe-type antiferromagnetic order below TV. Near the phase transitions, the system 
shows the s++-wave superconductivity due to the orbital fluctuation for a large g case with 
Tq > TV, while the s±-wave due to the magnetic fluctuation for a small g case with Tq < TV. 
The former case is consistent with the phase diagram of doped iron pnictides with T 3 > TV. 
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The recently discovered iron pnictide superconduc- 
tors 1 ^ RFePnOi-zF,, (R=Rare Earth, Pn=As, P) with 
a transition temperature T c exceeding 50K have at- 
tracted much attention. The parent compounds with 
x = show the tetragonal-orthorhombic structural tran- 
sition at T s and the stripe-type antiferromagnetic (AFM) 
transition at Tn . The carrier doping x suppresses both of 
the transition temperatures T s and Tn and induces the 
superconductivity. In HFePnOi- x F x , T s is always higher 
than Tjv, while in Ba(Fei_ a; Co a ;)2As2, the simultaneous 
first-order transition for nondoped case splits into two 



transitions with doping x where T s > T, 



N 



3) 



Theoretically, the s-wave pairing with sign change of 
the order parameter between the hole and electron Fermi 
surfaces (FSs), so called s±-wave, mediated by the AFM 
fluctuation was proposed as a possible pairing state in 
the iron pnictides. 4 ~ 6 ) The s±-wave state with a full su- 
perconducting gap seems to be consistent with most of 
the experiments. 7 ^ As for the impurity effects, however, 
the small T c -suppression against nonmagnetic impuri- 
ties 8 ' 9 ) is not consistent with the s±-wave state where 
T c is considered to rapidly decrease with the nonmag- 
netic impurities. 10 ) Therefore, the s-wave state without 
sign change of the order parameter, so called s ++ -wave, 
mediated by the orbital fluctuation which is enhanced 
due to the effects of the inter-orbital Coulomb interac- 
tion was proposed on the basis of the one-dimensional 
two-band Hubbard model 11, 12 -* and the two-dimensional 
16-band d-p model. 13 ) 

Remarkably, drastic softenings of the elastic con- 
stants have been observed in recent ultrasonic experi- 
ments. 14 ~ 16 ) As the elastic constant C e is given by the 
second derivative of the total energy w.r.t. the strain field 
e and includes the contribution such as —g^Xv with the 
susceptibility Xy for the electric operator f] linearly cou- 
pled with the strain field as g^jye, the enhancement of 
X v is responsible for the softening of C e . The detailed 
ultrasonic measurements 16 -* revealed that nondoped and 
underdoped Ba(Fei_ a; Co a: )2As2 shows drastic softenings 
of the elastic constants with decreasing T down to T s 
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Fig. 1. (Color online) Schematic figures of the strain fields for 
Ce, C44 and Ces modes (a), (b) and (c), the orbital fluctu- 
ations coupled with the corresponding strain fields (d), (e) and 
(f), and the phonons for B\ g , E g and orthorhombic modes which 
enhance the corresponding orbital fluctuations (g), (h) and (i), 
respectively. The x' , y' (x, y) axes are directed along the nearest 
(second nearest) Fe-Fe bonds. 



for Ce = (Cn — Ci2)/2, C44 and Cqq modes relevant to 
the strain fields e xx — £ yy , £ yz (szx) and e xy shown in 
Fig. 1 (a), (b) and (c) which are linearly coupled with 
the orbital fluctuations shown in Figs. 1 (d), (e) and (f), 
respectively, where x',y' (x,y) axes are directed along 
the nearest (second nearest) Fe-Fe bonds. The soften- 
ing of Cqq is much larger than C e and C44 and exhibits 
divergent behavior when approaching a critical temper- 
ature Tq which is just below T s . Then, the orbital sus- 
ceptibilities relevant to 6*66 mode, i. e., the longitudi- 
nal d y i z -d zx ' and/or transverse d 3z 2_ r 2-d xy shown in Fig. 
1(f), are considered to diverge at Tq below which the or- 
bital order occurs and induces orthorhombic distortion 
via electron-lattice coupling resulting in the structural 
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transition at T s ~ Tq. In addition, optimally doped 
BaFe1.84Coo.i6-A.s2 shows a significant softening of Cqq 
mode down to T c . 15 ' Therefore, the orbital order and its 
fluctuations relevant to Cqq mode are considered to play 
crucial roles in both the structural transition and the 
superconductivity. 

The orbital fluctuation is known to be enhanced by 
the electron-phonon interaction in addition to the inter- 
orbital Coulomb interaction. Recently, the effects of the 
electron-phonon interaction with B\ g and E g modes on 
the orbital fluctuation and its induced s++-wave super- 
conductivity have been investigated on the basis of the 
5-band Hubbard model 17 ' and the 16-band d-p model. 18 ' 
As shown in Fig. 1, the B\ g phonon enhances the lon- 
gitudinal d yz -d zx and transverse d 3z 2_ r 2-d x 2_ y 2 orbital 
fluctuations responsible for the softening of Ce mode, 
while the E g phonon enhances the transverse d x 2_ y 2-d yz , 
d X y-d zx and d 3z 2_ r 2-d yz orbital fluctuations responsible 
for the softening of C44 mode. However, the effect of 
the orthorhombic mode which enhances the longitudi- 
nal d y > z -d zx > and transverse d 3z 2_ r 2-d xy orbital fluctua- 
tions responsible for the most dominant softening of Cqg 
mode was not considered there. 17,18 ' The present paper 
is a straight-forward extension of our previous work 18 ' to 
include the orthorhombic mode which enable us to repro- 
duce the ultrasonic softening of Cqq and to obtain the x- 
T phase diagram including the tetragonal-orthorhombic 
structural transition and the superconductivity. 

Our Hamiltonian of the two-dimensional 16-band d- 
p Holstein model, in which 3d orbitals (d 3z 2_ r 2, d x 2_ y 2, 
d X y, d yz , d zx ) of two Fe atoms (Fe 1 =^4, Fe 2 =i3) and 4p 
orbitals (p x , p y , p z ) of two As atoms are explicitly in- 
cluded, is given by 18 ' 



(a) 



(b) 



(1) 



where H , H lnt , H p h and H e i_ p h are the kinetic, 
Coulomb interaction, phonon and electron-phonon inter- 
action parts of the Hamiltonian, respectively. The kinetic 
part of the Hamiltonian H n includes the atomic energies 
and the transfer integrals which are determined so as 
to fit both the energy and the weights of orbitals for 
each band obtained from the tight-binding approxima- 
tion to those from the density functional calculation for 
LaFeAsO and are listed in rcf. 13. In the present model, 
the doping concentration x corresponds to the number of 
electrons per unit cell n = 24 + 2x and there are two hole 
FSs (FS1 and FS2) around the T point and two electron 
FSs (FS3 and FS4) around the M point for x = 0.1. The 
Coulomb interaction part H m t includes the multi-orbital 
interaction on a Fe site: the intra- and inter-orbital di- 
rect terms U and V, Hund's rule coupling J and the 
pair transfer J'. For simplicity, we assume the relation 
U = U' + 2J and J = J' throughout this paper. Here- 
after, we number the Fe-3d orbitals as follows: d 3z 2_ r 2 (1), 
d x 2_ y 2(2), d xy (3), d yz (A), d zx (5). 

Now we consider the effect of the phonon and the 
electron-phonon interaction parts of the Hamiltonian 
Hph and H e i- p h which includes the phonon energy w s 
and the electron-phonon coupling constant g l s l between 
the orbital I and I' (see Fig. 2 (a)), respectively, where 
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Fig. 2. Diagrammatic representation of the electron-phonon cou- 
pling g s (a) and that of the bare vertex for the charge-orbital 
susceptibility C (b). 



s represents the phonon mode. In the present paper, we 
consider the B\ g , E g and orthorhombic modes as shown 
in Figs. 1 (g), (h) and (i). We note that the orthorhombic 
mode is not a normal coordinate but a general coordinate 
which is given by a linear combination of normal coor- 
dinates including both optical and acoustic modes. To 
avoid the difficulty with many phonon modes, we treat 
the orthorhombic mode as a local phonon similar to the 
B\g and E g modes as a simplest first step in includ- 
ing the orthorhombic mode. As following refs. 17 and 
18, we take the electron-phonon coupling into account as 
the atomic energy variance of the Fe-3cZ electrons. The 
resulting electron-phonon coupling matrix elements are 
given as follows: V^g^i = #| 5 i = 9%\ = ~^9e\ = 
9%\ = -9% = 9E S , 9 £ = -9%\ 9 = V3/2^ 2 l9 = 9B lg , 

— V / 3/25g 3 = 5q 5 = ge, g E s l — g l s l and for otherwise, 
where E y and E y correspond to the oscillation of the Fe 
atom along the x- and y-axis, respectively, and 9 denotes 
the orthorhombic mode (see Fig. 1). 

Within the RPA, 19 ' the spin susceptibility X s (l) an d 
the charge-orbital susceptibility x c {l) are given in the 
50 x 50 matrix representation as follows, 6, 13, 18 ' 

( %), 



X s (q) 

x c (i) 



X (0) (q)SY 



[i 



[i + x (0) («)C']- 1 x (0) (?) 



(2) 
(3) 



1 a,(8 



with the noninteracting susceptibility, [x^ ' (<?)] . \ 
-{T/N)^ k Gf° ei (k)G^f u (k + q), where a, p\=A,B) 
represent two Fe sites, £ represents Fe 3d orbitals, G{k) = 
\{ie n + fj,)l — Ho(k)]^ 1 is the noninteracting Fe-3d elec- 
tron Green's function in the 10x10 matrix representa- 
tion, /i is the chemical potential, Ho(k) is the kinetic 
part of the Hamiltonian with the momentum k given in 
cq. (1), fc = (k,ie n ), q = {q,iv m ) and s n = (2n + 1)itT 
and v m — 2rmrT are the fermionic and bosonic Matsub- 
ara frequencies, respectively. It is noted that when the 
largest eigenvalue A spin (A C _ G ) of x^^S (-x (0) (</)<?) 
reaches unity, the magnetic (charge-orbital) instability 
occurs. In eqs. (2) and (3), bare vertices for the spin and 
charge-orbital susceptibilities S and C are given by 17, 18 ' 

I Q\ot,fi _ (f T s\a.JS (r<\ a >P — lTTC\ a ^ 

2<W E s 9l 2ei 9 e s 3e4 D s (iv m ), where D s {iv m ) = 2u s /{v 2 m + 
uig) is the local phonon Green's function for the mode s 
(see Fig. 2 (b)). The bare vertices due to the Coulomb 
interaction tM c ' are given by, (U s{c) )"" u = U (U), 

(*> (C) )XV = U ' (-U' + 2J), {U'U)% t , t ,'= J (2U'-J) 
and (U s ( c) )%? e , e = J' (J'), where I ± t' and the other 
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Fig. 3. (Color online) Several components of the charge-orbital 
susceptibility x c («5f, 0) for U' = 1.0, J = 0.2 and g = 0.065 at 
x = 0.1 and T = 0.036, where we number the Fe-3d orbitals as 



follows: 



i(l), d a! 2_„2(2), <fe B (3), d yz (4), d zx (5). 



matrix elements are 0. 

The linearized Eliashberg equation is given by 

T 



AscA^(fc) 



N ^ 

k' l.xl.2 



k') 



r~iQL a ( 



a',/3' 



(k')G p J(k'), 



(4) 



where A(fc) is the gap function and V(q) is the effective 
pairing interaction for the spin-singlet state. Within the 
RPA, V(q) is given in the 50 x 50 matrix, 



V{q) = lsx s {q)S-\cr{q)C- 



5( A + ») 



(5) 



2 /v - v " 2 /v w 2 
The linearized Eliashberg equation (4) is solved to ob- 
tain the gap function A(fc) with the eigenvalue A sc . 
At T = T c , the largest eigenvalue A sc becomes unity. 
We use 32 x 32 k point meshes and 512 Matsubara 
frequencies (— 51l7rT < e n < 511ttT) in the numeri- 
cal calculations for eqs. (2)-(5). For simplicity, we set 
wb 1s = W£i = Wga = tug = ujq = 0.02eV as done in 
the previous study. To reproduce the experimental re- 
sults that the elastic softening is the largest for the Cqq 
mode, 16 ' we assume gB lg — gE g — 0-85gg and put gg = g. 
Here and hereafter, we measure the energy in units of eV. 

Fig. 3 shows several components of the static charge- 
orbital susceptibility x c (q,0) for U' — 1.0, J = 0.2 and 
g = 0.065 at x = 0.1 and T = 0.036. In this case, the di- 
mensionless electron-phonon coupling parameter is given 
by A = 2g 2 p /uJo ~ 2g 2 /ujq — 0.42 with the density of 
states at the Fermi level po ~ 1/eV. We hnd that, when 
T decreases, the transverse d yz -d zx orbital susceptibil- 
ity [x c (q, 0)]f 5 ' 4 5 , which is equivalent to the longitudinal 
dy'z-dzx' one, is most enhanced as compared to the other 
orbital and magnetic susceptibilities (not shown) due to 
the cooperative effects of the electron-phonon interac- 
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Fig. 4. (Color online) Several components of the gap function 
A(fc,i7rT) for U' = 1.0, J = 0.2 and g = 0.065 at x = 0.1 
and T = 0.036 (a)-(d), and those for V = 1.48, J = 0.2 and 
g = 0.032 at x = 0.1 and T = 0.034 (e)-(h). 



tion with the orthorhombic mode and the inter-orbital 
Coulomb interaction JJ'. 13 ' 1S ' We note that the incom- 
mensurate peaks around q = (0, 0) largely depend on the 
electron-phonon coupling strengths gj- and move to the 
commensurate peak at q — (0,0) for a slightly different 
parameter set where the resulting pairing state and the 
phase diagram discussed below is essentially unchanged. 

In Figs. 4(a)-(d), we show several components of 
the gap function with the lowest Matsubara frequency 
A(fc, ittT) obtained by solving the linearized Eliashberg 
equation (4) for the same parameters as in Fig. 3. In 
this case, the enhanced orbital susceptibility x c (l) f° r 
q ~ (0,0) (see Fig. 3), i. e., the ferro-orbital fluctuation 
yields the large negative value of the effective pairing 
interactions V(q) for q ~ (0,0) due to the 2nd term of 
r.h.s. in eq. (5), resulting in the gap function without sign 
change, i. e., the s++-wave state. For comparison, we also 
show the gap function in the case with a smaller (larger) 
value of g (U'), U' = 1.44, J = 0.2 and g = 0.032, in Figs. 
4(e)- (h). In this case, the enhanced magnetic susceptibil- 
ity x s (q) for q ~ (tt, it) (not shown), i. e., the stripe- type 
AFM fluctuation yields the large positive value of V(q) 
for q ~ (tt, tt) due to the 1st term of r.h.s. in eq. (5), re- 
sulting in the gap function with sign change, i. e., the s±- 
wave state. We note that, the effects of the ferro-orbital 
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Fig. 5. (Color online) Phase diagram on the x-T plane for U' = 
1.0, J = 0.2 and g = 0.065 (a) and that for U' = 1.48, J = 0.2 
and g = 0.032 (b). The symbols represent the instabilities for 
the ferro-orbital order (diamonds), the stripe- type AFM order 
(squares), the s+_|_-wave superconductivity (open circles) and the 
s±-wave superconductivity (closed circles), respectively. 



and the AFM fluctuations on the superconductivity do 
not compete to each other as they are mainly responsible 
for the different q regions in V(q), in contrast to the case 
with the antiferro-orbital and the AFM fluctuations. 

Fig. 5(a) shows the phase diagram on the x-T plane 
in the case with a large g, U' = 1.0, J = 0.2 and 
g = 0.065. When T decreases, the orbital susceptibili- 
ties {x c (q,Q)}fi% and [x c (.q,0)]?i% with Q ~ (0,0) (see 
Fig. 3(c)) diverge at a critical temperature Tq. Below Tq, 
the ferro-orbital order with different occupations of the 
d y i z and d zx ' orbitals occurs and induces the orthorhom- 
bic distortion resulting in the tetragonal-orthorhombic 
structural transition at T s = Tq . When approaching Tq , 
the ferro-orbital fluctuation is largely enhanced and me- 
diates the s ++ -wave superconductivity (see Figs. 4(a)- 
(d)). We also investigate the same model within the 
Hartree-Fock approximation 20 -* and obtain the phase di- 
agram consistent with RFePnOi-^F^ where T s — Tq is 
always higher than the stripe- type AFM transition at Tv 
for the case with relatively larger (smaller) orbital (mag- 
netic) fluctuation, and also obtain the phase diagram 
consistent with Ba(Fei_ :! ;CO; ! ;)2As2 where the simultane- 
ous first-order transition T s = Tq = T/v for x — splits 



into two transitions T s — Tq > T/v with doping x for 
a relatively smaller (larger) orbital (magnetic) fluctua- 
tion case. For the both cases, the ferro-orbital fluctuation 
dominates over the AFM fluctuation above T s — Tq for 
x > where the s ++ -wave superconductivity is realized. 

For comparison, we also show the x-T phase dia- 
gram in the case with a smaller (larger) value of g (f 7 ), 
U' = 1.44, J = 0.2 and g = 0.032, in Fig. 5(b). When T 
decreases, the magnetic susceptibility with q ~ (jr, it) di- 
verges at T/v below which the stripe-type AFM order oc- 
curs and induces the ferro-orbital order 21 ) together with 
the orthorhombic distortion resulting in the tetragonal- 
orthorhombic structural transition at T s = Tv, although 
the RPA result of Tq is smaller than that of Tv • When 
approaching T/v, the AFM fluctuation is largely en- 
hanced and mediates the s±-wave superconductivity. 4 ~ 6 ) 
In this case, the simultaneous phase transition takes 
place at T s = Tv even for x > and is inconsistent with 
the phase diagram of doped iron pnictides with T s > T/v 
which is reproduced for a large g case mentioned above. 

In summary, we have shown that the ferro-orbital fluc- 
tuation relevant to the ultrasonic softening of Cqq is en- 
hanced by the electron-phonon coupling g with the or- 
thorhombic mode and diverges at Tq below which the 
ferro-orbital order with different occupations of the d y i z 
and d zx i orbitals occurs and induces the orthorhom- 
bic distortion resulting in the tetragonal-orthorhombic 
structural transition at T s — Tq . Near the transition, the 
s ++ -wave superconductivity is realized due to the ferro- 
orbital fluctuation. The obtained x-T phase diagram is 
consistent with the phase diagram of doped iron pnic- 
tides with T s > T/v, in contrast to a relatively small g 
case with T s = Tv where the s±-wave superconductiv- 
ity is realized due to the antiferromagnetic fluctuation. 
For both cases with s++- and s±-wave superconductiv- 
ities, the RPA result of T c is always higher than that 
of Tq and/or T/v, where the orbital and/or the mag- 
netic fluctuations diverge. With including the effects of 
the self-energy correction and the vertex correction ne- 
glected in the RPA, it is expected that the ferro-orbital 
order and/or the antiferromagnetic orders are realized for 
relatively small x, while the superconductivity is realized 
for relatively large x. The explicit caluculation including 
such effects is a future problem. 

The obtained s++-wave superconductivity seems to be 
consistent with experimental results of iron pnictides in- 
cluding the impurity effects. The enhanced ferro-orbital 
fluctuation above Tq might be observed by experiments 
with a kind of external field inducing the anisotropy of x' , 
y' axes, similar to the case with the ferromagnetic fluctu- 
ation above the Curie temperature observed by experi- 
ments with the external magnetic field. In fact, a resistiv- 
ity anisotropy for T > T s is induced by uniaxial stress. 22 ) 
In the present paper, we treated the orthorhombic mode 
as a optical phonon, as a simplest first step. More real- 
istic model including acoustic phonons together with a 
suitable parameter set of the electron-phonon coupling 
strengths g% will be discussed in a subsequent paper. 
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